Dynamics of two interacting Bose-Einstein condensates 
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We analyze the dynamics of two trapped interacting Bose-Einstein condensates in the absence of 
thermal cloud and identify two regimes for the evolution: a regime of slow periodic oscillations and a 
regime of strong non-linear mixing leading to the damping of the relative motion of the condensates. 
We compare our predictions with an experiment recently performed at JILA. 
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The experimental evidence for Bose-Einstein conden- 
sation in trapped atomic gases Jl| has attracted a lot 
of attention, as the presence of a macroscopically occu- 
pied quantum state makes the behavior of these gases 
drastically different from that of ordinary gas samples. 
Trapped Bose-Einstein condensates are well isolated from 
the environment and, at the same time, can be excited 
by deforming the trap or changing the interparticle in- 
teraction. The question of how the gas sample, being 
initially a pure condensate, subsequently reaches a new 
equilibrium state is directly related to the fundamental 
problem of the appearance of irreversibility in a quan- 
tum system with a large number of particles. Thus far 
the time dependent dynamics of trapped condensates has 
mainly been analyzed for a single condensate |^^| on the 
basis of the Gross-Pitaevskii equation for the condensate 
wavefunction. Remarkably, already in this mean field 
approach the stochastization in the condensate evolution 
has been found ||, and the damping of the condensate 
oscillations has been observed numerically 0. However, 
the question of the formation of a thermal component, 
addressed in has not been investigated. 

In this paper we study the evolution of a richer sys- 
tem, a mixture of two interacting condensates (a and b) , 
in the situation where initially the thermal cloud is ab- 
sent. The properties of a static two-component trapped 
condensate, including the issue of spatial separation of 
the a and b components due to interparticle interaction 
M,|J, were investigated in ||. The response of the sys- 
tem to small modulations of the trap frequency has also 
been studied numerically |lC|| . In our case the a and b 
condensates have initially the same density profile and 
are set into motion mostly by an abrupt displacement of 
the trap centers. The main goal of our work is to study 
the dynamics of spatial separation of the two conden- 
sates and analyze how the system can acquire statistical 
properties and reach a new equilibrium state. From a 
general point of view, we are facing the problem raised 
by Fermi, Pasta and Ulam [jnj. They considered classical 
vibrations of a chain of coupled non-linear oscillators, to 
analyze the emergence of statistical properties in a sys- 
tem with a large number of degrees of freedom. As has 
been revealed later, the appearance of statistical prop- 



erties requires a sufficiently strong non-linearity leading 
to stochastization of motion |l2j , whereas for small non- 
linearity the motion remains quasiperiodic (see e.g. Jul). 

We consider a situation in which the two condensates 
a and b see harmonic trapping potentials of exactly the 
same shape, and the interparticle interactions character- 
ized by the scattering lengths a aa , a a b and are close to 
each other. The control parameter, determining the pos- 
sibilities of non-linear mixing and stochastization, is the 
relative displacement z of the trap centers. We identify 
two regimes for the evolution. In the first one the rela- 
tive motion of the condensates exhibits oscillations at a 
frequency much lower than the trap frequency lu. In the 
other regime there is a strong non-linear mixing leading 
to the damping of the relative motion, and the system has 
a tendency to approach a new equilibrium state. We com- 
pare our predictions with the results of the JILA exper- 
iments [ |14|Jl^1 on a two-component condensate of 87 Rb 
atoms in the F = l,m = — 1 and F = 2,m = 1 states. In 
these experiments the double condensate was prepared 
from a single condensate in the state F = l,m = — 1 
(a) by driving a two-photon transition which coherently 
transfers half of the atoms to the state F = 2, m = 1 (b). 

We mostly perform our analysis in the mean field ap- 
proach relying on the Gross-Pitaevskii equations for the 
wavefunctions <f> a and fa of the a and b condensates. This 
approach corresponds to the classical limit of the evolu- 
tion of a quantum field, the subsequent corrections be- 
ing proportional to a small parameter (na^,) 1 / 2 (n is the 
gas density) and, hence, manifesting themselves only on a 
rather large time scale. The two coupled Gross-Pitaevskii 
equations for 4> a and fa normalized to unity read 
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Here g ee i = Airh a ££ i /m are the coupling constants for 
elastic interaction between atoms in the states e and 
e', m is the atom mass, and N ei U £ are the number of 
atoms and trapping potential for the e condensate. As 
in the JILA experiment, we choose the initial condition 
fa,b(0) = fa, where the (real) wavefunction fa corre- 
sponds to the ground state of Eq.(|]) with all atoms in 
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the a state and no trap displacement. The chemical po- 
tential of this ground state is denoted as /i. 

We consider the a and b condensates in the Thomas- 
Fermi regime (Kuj <C and assume the number of con- 
densate atoms N a =N b = N/2 (r|. The first set of our 
calculations is performed for the evolution of the con- 
densates in a spherically symmetric trapping potential 
Uo(r) = muj 2 r 2 /2 which at t = is displaced along the z 
axis by a distance zq/2 for the a atoms, and by —zq/2 for 
the b atoms. We present the results for the time depen- 
dence of the mean separation between the condensates, 
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For the curves in Fig.l the coupling constants are g aa = 
9ab = 9bb, and for zq = our initial state is an equilib- 
rium state at t > 0. In this state the Thomas-Fermi ra- 
dius of the condensate Rq = {2^,/muj 2 ) 1 / 2 serves as unit 
of length, and the shape of 4>q is determined by [i/Tiuj. 
Hence, for z ^ the dependence of the quantity u/Rq 
on wt is governed by the parameters ji/hui and Zq/Rq. 




FIG. 1. Mean separation between the condensates ver- 
sus time in isotropic traps for g aa = g a b = gtb and 
/i/hui — 29.2. Relative displacement: zq = 6.66 x 10 _4 7?o 
(a), and zo = 7.16 x 10" 2 i?o (b). Solid curves: numerical in- 
tegration of Eq.(|l]). Dotted curves: analytical prediction for 
(a) (see text), and the linear model relying on Eq.(^|) for (b). 

Our results reveal two key features of the evolution 
dynamics. The first one, for a tiny displacement zq, is 
a periodic motion with slow frequencies which turn out 
to be sensitive to small variations in the values of the 
coupling constants. The other feature, for much larger 
Zq, but still zq <C Ro, is a strong damping in the relative 
motion of the two condensates, as observed at JILA p4[ . 

In order to understand the physics behind the evolu- 
tion pattern, we first perform a linear analysis of Eq. ([l]) . 
For the case where g aa = g a \, = g^ = <?, and the displace- 
ment zq is sufficiently small, we linearize Eq.(fJ) with re- 
spect to small quantities 5<j) a ,b — {4'a,b—4 l o) an d zq. Then, 
for the quantity 5<fi- —8(j) a — 84)b, describing the relative 
motion of the condensates, we obtain the equation 



H 2 A 
2m 



+ Uq-li + Ngc/) 2 



(3) 



with the source term S- = rmo 2 ZqZ<\>q. For the quantity 
5tf>+ = 6<fi a + 8<fib we find an equation decoupled from 54>- 
and without source terms. Hence, the initial condition 
5<j)+ (r, 0) = allows us to put 5<fi+(r, t) — for t > 0. 

For S- = Eq. (||) is the equation for the wavefunction 
of a particle moving in the potential V = Uq — n+Ng^. 
Stationary solutions of this equation provide us with the 
eigenmodes of oscillations of the condensates with respect 
to each other. In the Thomas-Fermi limit the potential 
V, originating from the kinetic energy of the condensate, 
is a smooth function of r inside the condensate spatial re- 
gion r < R Q : V = h 2 (A(/)Q)/2m(f)Q < hto. For r > R this 
potential is close to Uq — fi and is much steeper. Replac- 
ing V by an infinite square well of radius Rq we obtain 
the energy spectrum of eigenmodes with large quantum 
numbers n: E n j = (itTujj) 2 {2n + l) 2 /16fi, where I is the 
orbital angular momentum. This explains the appear- 
ance of oscillations at a frequency much smaller than u 
in our numerical calculations (see Fig. la), since the en- 
ergy scale in the spectrum is (Ucu) 2 / ' \l <C Uco, For the 
latter reason we call these eigenmodes soft modes. Note 
that the soft modes for the relative motion of the two 
condensates also exist in the spatially homogeneous case 
where they have a free-particle spectrum 0. 

As in our linear approach we have 5(j)+(v, t) = 0, Eq.(||) 
for the mean separation between the condensates reduces 
to u(t) — 2 J d 3 r z (f)Q ~Re{6(j)-}, and the contribution to 
u(t) comes from the components of S(f>- with I — 1, m; = 
0. Solving Eq.(|J) with the initial condition <50_ (r, 0) = 0, 
we obtain u(t) as a superposition of components, each of 
them oscillating at an cigcnfrcquency of a soft mode: 
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where ip n \ is the wavefunction of the soft mode with 
I = 1, mi = and main quantum number n. Damping of 
oscillations of u(t) could, in principle, originate from the 
interference between the components with different n in 
Eq.(|j). However, the source S- basically populates only 
the lowest soft mode, irrespective of the value of zq: the 
amplitude of oscillations at the lowest eigenfrequency in 
Eq.(Q) (the term with n = 1) greatly exceeds the sum of 
the amplitudes of other terms. Hence, these oscillations 
remain undamped. For the same reason their frequency 
and amplitude can be found with tp n i replaced by the 
function z(f>Q normalized to unity. Using the Thomas- 
Fermi approximation for the condensate wavefunction 
0: $(r) = 15(1 ~r 2 /R 2 )/8nRl for r<Ro, and O = 
for t>Rq, we obtain En = hil — (7 / 'A)(fiui) 2 / ' \i which is 
very close to E\\ = \ .62{hio) 2 /n calculated numerically. 
Then, retaining only the leading term (n = 1) in Eq.(||), 
we find u(i) « zo(4/i/7?iijj) 2 [l — cos(51t)] shown in dotted 
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line in Fig. la. As one can see, the condition of the linear 
regime u <C Rq requires a very small displacement 



z < (huj/ii) 2 R , 



(5) 



and already a moderate zq as in Fig. lb is sufficient to 
drive the system out of the linear regime. 

We have performed a similar linear analysis for the 
case where g aa g a b ^ 9bb, but the relative difference 
between the coupling constants is small. Also in this case 
the source SL mostly generates oscillations of the conden- 
sates relative to each other at a single frequency SI' <C ui. 
For a relative difference between the coupling constants 
much smaller than (hcu/ii) 2 , the frequency SI' coincides 
with the soft-mode frequency SI found above. Otherwise 
the sign of </_ = g aa + gbb — 2g a b becomes important. In 
particular, for positive g_ ^ \g aa — gbb\ already a mod- 
erate difference between the coupling constants strongly 
increases the frequency 51' compared to SI. In this case 
we obtain undamped oscillations at SI' ss (g- / ' g aa ) x ^ 2 u>. 
For g_ < 0, already in the z = breathing 
mode in which the two condensates oscillate out of phase 
becomes unstable, and the system evolves far from the 
initial state. Note that for a small difference between the 
coupling constants the condition <?_ < is equivalent to 
the criterion of spatial separation of the condensates in 
the homogeneous case, g aa gbb < g 2 ab @§] ■ 

We now turn to the large zq regime (Fig. lb) where 
we find a strong damping of the oscillations of the mean 
separation between the condensates, u(t). In order to 
prove the key role of non-linearity in this regime, we first 
attempt a linear model assuming that the densities \<p £ i \ 2 
inside the square brackets of Eq. (fy) are not evolving: 
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In contrast to the analysis which led to Eq. ([!]), the dis- 
placement zq is now explicitly included in the Hamil- 
tonian through the terms ±muj 2 zzq/2 in U a ,b, and the 
number of populated oscillation modes depends on z . 
However, for the parameters in Fig. lb we find that only 
a few modes are populated, and the interference between 
them can not account for the damping found numerically 
(dotted versus solid curve in Fig. lb). 

We argue that the damping in our calculations mostly 
originates from non-linearity of the system, which in- 
creases the number and amplitude of populated oscilla- 
tion modes and provides an interaction between them. 
As a result, the evolution of the condensate wavefunc- 
tions 4>a and tpb becomes chaotic. This can be seen from 
Fig. 2 where we compare the spectral density i?„(^) — 

IT' 1 J Q dtn(0,t) exp(ii't)\ 2 of the density at the origin 
n(0, t) with an identically defined spectral density R u {v) 
of u(t) for the parameters in Fig. lb and T = 110/w. 
The function R n {v) has a smooth envelope at large 
with peaks corresponding to the islands of regular mo- 
tion. On the contrary, R u (v) exhibits pronounced peaks 



at v of order u>, without any smooth background. This 
picture provides a clear signature of stochastization in 
the system Q and prompts us to represent each of the 
condensate wavefunctions in Eq.(Q) as a superposition 
of two constituents: (i) a slowly oscillating regular part 
conserving the phase coherence properties; (ii) a composi- 
tion of high-energy excitations characterized by stochas- 
tic motion. Only the slow constituent contributes to such 
macroscopic quantities as u(t), since the contribution of 
the fast stochastic part is averaged out. 
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FIG. 2. Spectral densities R„(v) (a) and Ru(v) (b) for 
the parameters in Fig. lb and T = 110/w (see text). 

Our analysis is consistent with the general statement 
that for a large population of various oscillation modes 
the non-linear interaction between them leads to stochas- 
tization in the motion of excitations with sufficiently high 
energy This allows us to employ the mechanism of 
stochastic heating |H| for explaining the damping of os- 
cillations of u{t): The mean field interaction between the 
fast stochastic and the slowly oscillating parts leads to 
energy transfer from the slow to the fast part. 

The evolution of the occupation numbers of the modes 
of the fast stochastic part is governed by kinetic equa- 
tions Jl3| ] and eventually slows down. The rate of energy 
and particle exchange between the two constituents then 
reduces. After a sufficiently long time only small linear 
oscillations of the condensates survive, mostly at the low- 
est eigenfrequency and the gas sample as a whole could 
be thought as being close to a steady state. However the 
damping of the remaining oscillations and the ultimate 
evolution of the fast stochastic part towards the thermal 
equilibrium require an analysis beyond the mean field 
approach. For the parameters in Fig. lb, using the semi- 
classical Bogolyubov approach Um and relying on the 
conservation of energy and number of particles, we find 
an equilibrium temperature T oq ks 0.6/i and a condensed 
fraction la , b w 0.9, for TV = 5 x 10 5 . 

The last set of our calculations relates to the recent 
JILA experiment juj where the evolution of a two- 
component 87 Rb condensate has been investigated. In 
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the conditions of this experiment we solved numerically 
Eq.(||) by taking a ab = 55A and the ratio g aa : g ab : g bb = 
1.03 : 1 : 0.97. We also explicitly included in these equa- 
tions the 22 ms expansion of the clouds after switching 
off the trapping potential. The results of our calculations 
are presented in Fig. 3. As in Fig. lb, we find a strong 
damping of the oscillations of the mean separation be- 
tween the condensates, u(t). Our numerical results are 
in fair agreement with the experimental data, although 
the damping in the experiment is somewhat larger. We 
extended the calculations to twice the maximum exper- 
imental time and found small oscillations which remain 
undamped on this time scale. 
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FIG. 3. Mean separation between the condensates in the 
JILA experiment versus evolution time in the traps, after a 22 
ms free expansion. Dots with error bars: JILA experiment. 
Solid curve: our numerical calculation. 

Our data for the JILA experiment GJJ] can be ana- 
lyzed along the same lines as the results in Fig. lb, with 
a damping originating from stochastization in the evolu- 
tion of the condensate wavefunctions. The equilibrium 
temperature is close to //, corresponding to condensed 
fractions 7 a rs ^ b w 0.9. The large value of the con- 
densed fraction explains why phase coherence between 
the a and b components could be observed even after the 
damping of the motion of u(t) pJJ. The damping time 
of the small remaining oscillations, estimated along the 
lines of p9| , will be of order 1 second. 

We believe that the stochastic regime identified from 
our calculations is promising for investigating the loss of 
phase coherence and the formation of a new thermal com- 
ponent in initially purely Bose-condensed gas samples. 
An interesting possibility concerns the observation of a 
continuous change in the phase coherence between the a 
and b components with increasing the trap displacement 
and, hence, decreasing the final Bose-condensed fraction. 
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